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ABSTRACT 



Tlie determination, af eTasiTc: properties of : crystals :by 
sound velocity measurements^ is: f^iiitated by use-of ;ai 
perturbation soIutTon. to the; e.quations of motion: for; sound 
waves propagating in. the cnystaiLs:. . Transformation; equations 
can relate the solution fbcr- arbitrary porpagat ion . direction 
to an established crystai eojordlhate system'. A.-. method :of 
calculating the acoustic components- of the energy flow vector 
Is developed for severai high symmetry crystallclasses . The 
validity of this method', is- hereiir. confirmed and. results are 
presented which, are consistent; with- experimental '.data . . A 
modern desk; calculator' vjas employed; to determine-, solutions 
for hexagonal,, cubic and; tatragonail crystals. . 
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A major portion of the fehecEretieal dev eiopTnent. presented 
in this paper has previausly.' heen. presented:' by Priofi . John R. 
Neighbours, Physics Department,, U.. S.. Nav all Post graduate? 
School. The most recent paper,. ao:-authored:; by Associate? 
Professor G. E. Schacher,. Ehysi'ca D:epartment',, D.. . S.. . Naval 1 
Postgraduate School, is titled' ''Deterrainatxdn'. of; Elastic : 
Constants from Sound Yelocity Measurements' in: Cry.stal's; of : 
General Symmetry.”' In addition.,, several, investigators; have 
published methods similar to this development... These; papers 
are included in the List of Ref'erences^. 
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I. INTRODUCTIOIT 



The task of relating velocity nreasuireraentis t.a elars±ijC- 
constants of a crystal is complicated by tire fact that: tiie 
eigenvectors which describe the direction, of vibration in 
response to the sound wave are neitiier perpendicular nor. 
parallel to the propagation direction.. The energy flow may, . 
however, be described by modes of vibration, which are: quasl--- 
longitudinal (essentially parallel to the propagation: direc- 
tion). In practice one vibrational mode dominates the.; others: 
for a given propagation direction.. In addition, the: transducer: 
used to measure the velocity has greater response to: one: mode: 
of vibration so it is possible to separate modes by appropriater 
choice of transducers. 

For certain high symmetry crystal classes specla:! dlr.ec^- 
tions of propagation exist in which the response 'is pure 
longitudinal and pure transverse. Telocity measurements^ along; 
these directions can be directly related to cry stall elastic, 
constants. However the energy flow is not necessarily parallel' 
to the propagation. 

For either high or low sjnranetry cases the elastic: constants: 
can be determined by solution of the crystal vibration equation 
of motion. This process may be simplified by appropriate 
choice of coordinate system and application of several assump- 
tions. The resultant secular equation is solved by a 
perturbation method and the values for the elastic: constants: 
may be refined by Iteration. The constants thus obtained are 
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linear camhxnatlana clC the erl'aslr±c: constants; defined _by ; 
conrenti-onal c:acard:ina-t.e- syafceins: and: carr. bee converted :by ; 
rotational tacansf"cn2matxcaT^, The: angiee of: dev±ation::of I'the 
energy flow vector from the: propagation direction, is :of : 
particular interest in. order to. constrxict' a: physical Ipicture 
of the crystal response:.. This; angle; can be-' obtained :by 
application of the general energy flow equations : to 'the.-, 
coordinate system defined: by the,- propagation, direction. . 
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II. SECULAR EQUATION 



The general form of Hooke's Law Is 



( 1 ) 



Two common conventions of subscript notation for the above 
equation are shown below. 




The strain, S, is related to displacement by the equation 



\l 2 ^ ^k,l ^l,k^ 



1 f 

2 ^ au^ 



3U. 

aul 



) ( 2 ) 



where u is displacment. In contracted notation Hooke’s Law 
may be written as 



'^i “ ^ij^J 



( 3 ) 



where the following notation applied: 



= Li = 
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y X 12 21 6; 



The existence of a strain energy function, re-duces.'' the;- 3.6 
elastic constants to 21 Independent, constants: for.' the: most 
general case. This number can be further.' reduced: fnr.' specific 
crystal calsses because of symmetry considerations:. . 

If an arbitrary coordinate system is- chosen, such, that . the 
x' axis lies along the direction of propagation the: equations 
of motion can be written in terms of the elasti'c constants in 
the arbitrary coordinate system and dispiacement. . In. this: 
development u’ , v'' , w'' are displacements' aXong: the: x' , . y;* , ■ 
and z* axes respectively^ The primed no.tation emphasizes the 
fact that choice of coordinate system is- arbitrary,. . The^ 
equations of motion are 



P^, = ma 



p — dx^dy'dz*’ 

at 




3T^ aT.g 3T^ 




^ 3y^ ^ 3z' 



3Tg 3T^ 



( 5 ) 
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If plane wave prQpagatxan. 1 st assiuiredi ttier partiallderlva- 



tives with respect ta y'' and. z:*' vanish', and: thee equation 
becomes 



^11 



3 U' 



3X 



,•2 



- + 









3x 



4_ P'l ■ 

2 ^ ^15; 



3 w * 
3x;' 



2J" 



( 6 ) 



Similarly, the equations. o£‘ mollon'. fhr/ the: remaining^directions 



are 



2;..- 



3 V 



2 : 



p _ f,,. iSll-t-pf 



+-01,.^ 






=" p 



3^w*' 

3tt^ 



=■ CTl', 






t- Cl^ 



2;.,' 



3 V 






2 ■ 

3^W' 



15 3 -^;, 21 " 55 : 3 ^, 2 . " 55 : 3 ^', 2 



An additional assumption, is: that: tha; solutions: are .periodic 
of the form 



u" 


< 

li 


COS' 


(kx'-tol) 




v" 


II 

> 


CCTS^ 


(kx^-tjt:) 


(7) 


w’' 


= A3 


COS 


( kx ' wt ) 





The assumed periodic solutions' are inserted in the equations 
of motion and the partial derivatives, are evaluated. The 
resultant equations may be further simplified by the relation 
V = u/k where V represents velocity .and is not to be confused 
with v', the displacement in the y' direction. The equations 
may now be written 
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0 ( 8 ) 
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The secular equatrlon. has: three-: reaX roots corresponding to 
one quasi-longitudinal and twa transverse modes: of: vibration. 
If the solutionis are assumed: ta he essentially; longitudinal 
or transverse with, small deviations from these directions 



the secular equation, yields: a- perturbation solution: : 



r ' 
^11 



= P 









2 

A A fT T 









Vil-=66) 



Note that the amplitudes of'the assumed periodic: solutions do 
not affect the calculation.. The remaining solutions 'to the 
perturbation equations are. 



C 



66 




pVt - 



^16' 



•'56 



^66 - 4 l - ^55 



pV^ - 



^15 



^56 



( 9 ) 



- "il - "66 



In particular, if the terms in square brackets vanish the 
solutions correspond to pure mode response and the elastic 
constants are determined directly from the measured velocities. 
When the square bracketed terms' are present velocity measure- 
ments along several direations yield an adequate number of 
equations to determine all six elastic constants of the secular 
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equation to flirst order .. Uiesee vadites; aree thenr. inserted in 
the square bracketed terms and the; soiiLtlbn" is: ref ined . The 
primed elastic constants are: linear combihatlonsrof I'those 
defined by the conventicmaOl coordihaire:- system. The 'rotational 
transformation relating- the two systems: is: developed : in 
Section IV. 

The amplitudes of" the- periodic: soiutibns ■ may be -determined 
from the equations of matxan.. IT.' the: ampditude: of :'one mode 
of interest is fixed at an. arbitrary value- the: remaining two 
amplitudes can be expressed in. terns: of: the elastic .constants 
and the measured velocities.. For.- exampde , . if:' the amplitude 
is set equal to unity the' longitudihall raoddr solution is 



A 



2 



A 



3 




t* 





45 - 



P-V‘: - C • 5 





4— 


=56- 1 


1 


=15 


1 


^4 -“55. 


1 

J 


pvf - 




=15 ’ 




“is 


1 


“56 _ 




PV2 - 0^5 


J. 


“l5 




- “56 


1 


“i6 . 





( 10 ) 
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irr.. eneroy; flow 



The components of the Energy Elow vector/ are - 
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E , = + r'X' + w’X^ 

x’ X y z 



E , = + r’T'' w’Y^ 

X y z 



( 11 ) 



E , = u*’Z*’ + r'Y*' +- w^X*' 
z’ X y z 



The X* component of this ve<rtor can be calculated/as; follows 



Ex’ 



ti •' C c ’’ * +- (T' ■ ^ ^ * + n * ) 

^ ^11 ax'' ^i6l 9.x.' ^15. 8x' ^ 



v*- c CIV +- cyr + c’ . 

I.D ax' 6.6. 9x' 5.6 9x.' 

4 ^r. ( rv r*- + C’ 



1 



(12) 



Note that the plane wa.ve assxunption has been applied. to yield 
only those terms which, contain, partial derivatiyes: with respect 
to X*. The transverse components* of the Energy Flow vector are 



"y' 



= u’ f C ' ' + c' * + C' ' ■) 

^ ^l6 9x' ^66 9x' ^56 9x' ^ 



+ V* fC ^ + C ^ + C* ^ ^ 

^ "" ^^12 ax' ^ ^26 ax' ^ ^25 9x' ^ 



'26 ax 



(13) 



+ JrT (C^ iiLl + r* ^ + r* ^ 

^ ax' ^ ^46 ax' ^ ^45 ax' ^ 



■‘"A. E. H. Love,. A Treatise on the Mathematical' Theory of 
Elasticity (New York, 1944), ch. 7. ~~ 
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On substitution of the periodic solutrons- the x;*' caomponent. 
becomes 

E^, = (uA^sin n [Cj^^C-kA^sin n) + CT^^X-hA^siir. it) 

+ C^^C-kA^sin n)] + oiA^ai'n n i~ siii- nO 

+ C^g(-kA2Sin n) + C^gC-kA^slh n)l 



where n = kx' - ut . 



The above form may be further slmpIiTled hy d'efihihg: 



The final form of the Energy Flov: components is shown below. 



+ (jA^sin n [ 



P 








13 



'’z- ' 4s''? ^ * °sA * '46 * '4s ^ 43>44 

(Cilj + C^g)A2A3 (15) 

The resultant Energy Flow vector determines deviation from 
direction of propagation as shown in the diagram below. Note 
that the angle 6 Is the compliment of the polarization angle. 




FIGURE 2 



lil 



IV. TRANSFORMATION OF ELASTIC CONSTANTS 



The desired elastic constants are defined by a fixed, 

unprlmed coordinate system which is related to the crystal 

basis vectors a, b, c in accordance with convention adopted 

2 

by the I.R.E. standards committee. Specific examples of this 
convention are provided in Sections V, VI and VII. The primed 
coordinate system is chosen so that the x’ axis lies in the 
x-y plane. This leads to a relationship between the two 
systems in terms of the direction cosines 

£ = sin 6 cos i}i 

m = sin 6 sin 4) (l6) 

n = cos 6 




FIGURE 3 



2 

standard Committee, ’’Standards on Piezoelectric Crystals,” 
Proc. I.R.E. 37 (September, 19 ^ 9 ), p. 1378 . 
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Tiie trann forma t.lan. edlonenta are gl.vea’-.ih:'. tabular: form 



below. 





X. 


y 


z: 


x;*' 


1 


m 


rr. 


y'- 


-1 

-ma 


la“^ 


o; 


a*' 


! . -1 
1 — Inoc 
i 


-1 

-mna 


a- 



(17) 



The elements have been reduc-eh. by use of: the: definition 



2 2 2 
a = + m . 



For examp le :: 



« f 1 « “i- 

= lx;* — raa-. y’ -ina. 5 



( 18 ) 



The transformation afT eXaetic constantsrmay be : obtained 
from Hooke *'s Law;; 



Ij, = CySj 



(19) 



If the transf ormationa of stness and straih'.are T! = 
and S* = the elast.la' (Constants are relatediby the 

expression 

( 20 ) 

where = g” and g la the matrix shown below. 



£2 


2 

m 


2 

n 


2mn 


2ln 


2lm 


2 -2 
m a 


»2 -2 
1 a 


0 


0 


0 


-2lma~^ 


2 2-2 
1 n a 


2 2-2 
mna 


Z 

a 


-2mn 


-2ln 


„ 2 -2 
Imn a 


-2 

imna 


„ -2 
-imna 


0 


1 


-m- 


/o2 2y -2 

-n(l -m )a 




-m^na~^ 


na 


m(X— 2n^)a ^ 


l(l:r2n^) 


-1- on -1 

-2 Imna 


-ima”^ 


ima 


0 


ina””^ 


„,ry -1- 

-nina 


/«2 2\ -1 
(1 -m )a 



The Energy Flow equations and the rotational. transformation 
may novx be applied to specific crystal classes. 
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Y. CRYSTALS OF HEXAG'ONAF SYMMETRY. 



In application to aerta^in. higlr. ay iranetry; crystal! classes : 
the method is simplified due t:a symmetry considerations :and : 
the fact that certain primed elasiic: c:onst ants; may; vanish. 

This leads to a partially dia gonad! aerular.' equation: and: the i 
existence of semi-pure modes of' vihra.t'ion. . For: hexagonal . 
crystals the I.R...E.. convention, defines: the: unprlmed: coordinate 
system such that the z-axis is parali'el! tn: the: crystal !basis : 
vector c, i.e. the sixfold symmetry axis:,, and: the: Xf-axiS : is : 
parallel to a» the two f aid rotation' axiSs . The: y;-axis ;iS : 
perpendicular to r and z ta preserve a right: handedicoordiriate -■ 
system. With this definition, the primed: elastic: constants: 
become 

C|i = t 2n.^a^(C^2: 44 ^ 

C»3 = n^a^C^ - ~ ^ 33 ^ +- (i:-2h;^)^C^H ^ 

=66 ' I “"(=11 - ^ 

= -na'^C^^ + na(I-2n )(C,'^- F2C^^) +-n^ctC 33 

=16 = =56 ' = 



Since and C^g are zero the existence of a semi-pure mode 
is implied. Physically this means that” the. vibrational '.response 
of the crystal lies along one of the primed ‘ system-' axes . for : 
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one mode (pure transverse) and the remaining modes are quasl- 
longitudlnal and quasi-transverse . The relative orientation 
of the eigenvector components to the propagation is shovm 
below. 




FIGURE 4 



Note that If were to vanish the secular equation 

would become diagonal and the response would be pure longi- 
tudinal or transverse. This situation may occur for materials 

possessing hexagonal symmetry. The existence of this acciden- 

q 

tal pure mode was predicted by Borgnis-^ for the case of 
propagation in a direction that makes an angle 0 with the c 
axis defined by 



tan 



2 



6 



^33 ~ ^13 ~ 
^11 " ^13 " 



(23) 



The above requirement corresponds to the fact that vanishes 

and the secular equation is diagonal in the perturbation 
solution. 

Equations (22) are sufficient to calculate the velocities 
of each of the modes from known values of conventional elastic 



q 

■^P. E. Borgnls, "Specific Directions of Longitudinal Wave 
Propagation in Anisotropic Media," Phys. Rev. 98 (August, 
1955), p. 1000. 
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constants. Tlie elastic constLarrts: ccfT mniier.ous; hexagpnall. 
crystals have already been, determined: and: the: application of. 
the perturbation method may be demans-trafred: by comparison of 
calculated mode velocities with measured: valiiesi . The: solu-- 
tion for energy flow, however,, requires nine: additlonall 

constants i.e. ^i3>- ^14 . = 

Some of these may vanish as a result, oi symmetry.. For: 
hexagonal symmetry and C."’ '2^^. are- alll identically 

zero. The remaining primed constants are 



42 ' “^42 '^^43 

= Cl - C^3'. - 



^25 ria(C ^2 — ^2^2^ 



(2ft) 



= -na(n^C^2. 2X11-2 h^')C^ 4 ) 



^ft6 " " 2 ~ ^12 ~ ^^ 44 -^ 



In addition, the fact that 46 and Cl^ are- zero: leads: to 
simplified forms of the amplitude^ equations.. With" these 
simplifications the equations for energy flow in hexagonal; 
crystals are 



4- ' 4i'‘i ^ 4e4 4 s''3 ^ 



Py, = Cc^g + C^ 2 )Aj^A 2 + (C ^5 + C 4 g)A 2 A 3 ( 25 ) 

4- = 45*1 + 45*1 + <4? ^ 6r3)*i.*3: ^U6^2: 



19 



Some of the terms in. the- ahcrv.e- equations; wl3.ll vanlsh.-^ when 
certain, roots of the seculaar equat±on~. are:- applied since- the 
amplitudes of the aasxxmedl perijodle response: VcLnish.; . In the 
notation used below the second subscript denotes' thef root ‘ of 
the secular equation that is used, iir the: calculation. . For 
instance, is the amplitude when, the' root: is.- used. 

The amplitudes are 

A=1 A— CTA-— — — ^ 2 ^ 

^11 ■ ^ ^12 ^13- C' - 

IP- 




A 



22 



£ L 









11 



( 26 ) 



The energy flow relationships' for.' the modes, of vibrational 
response in hexagonal materials are 

For the root (Quasd-CongltudinalO 



,.2 



».2 






P , = 0 

y 

^ ^ 15^11 ^ 35^31 ^^ 5 . 5 ; ‘^ 13 :^^ 11 - 31 - 



( 27 ) 



For the root (Pure Transverse) 



^x^ ^66^22 



Fy. =0 



F^,.- = C/,^A 



( 28 ) 



46 22 
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( 29 ) 



Par the raat ( Quasjl-Tr.ans3^er:ser) 

2 2 

^x:’’ " ^ix^L3 ^ ^^15^13^^33; 

= Q 

^z’- "" ^15^13"^ *^fe^33; ^13?^13^33 

For this symmetry the ahove- three equetihns: (27j -28; .and 29) 

show that the y*' campanent o£ the: energy flow vector, is 

always zero. Thus the ang^Ie tp. deiiiied* ort. page 1^- is. ‘always 

zero and any deviation o£' the- energy flow vector: from the 

propagation vector tahes pla.ee in the: r’-z:’' plane.- . 

In order to demonstrate the? validit'y.' of:‘ the-:-perturbation 

method solutions were computed.; on a? Hewlett*. Packard: 98 IOA 

desk calculator and compared- to: published: data: for. hexagonal 

crystals, specifically zihe and' beryl,. as:presented:by M. J. 

4 . 

P. Musgrave« Transformation equations: using, known. '.values of 
elastic constants were programmed and the;- velocities obtained 
were in agreement to ..Q5S.. The* energy • flow components were 
calculated and the angles, of deviation fronr propagation 
direction were consistent to within 2°' of. those presented by 
Musgrave. The graphs on the foU'owing pages are examples of 
results obtained for x-z plane propagation in hexagonal 
crystals.. In these graphs the deviation angle, 6 ,. Is plotted 
as a function of the angle,. 6 ,. which the propagation vector 
makes with the z-axis for each of the; three modes. In addition 



2i 

M. J. P.. Musgrave, ”0n the Propagation of Elastic Waves 
in Aeolotropic Media," Proc. Roy. Soc . 226 (195^) p.-356. 
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the sign of all deuxatiorr. angles' in the:- cubic: and:l hexagonal 
cases has been reversed in: ond'en' ta ceinpane: the-: results with 
Musgrave's data* Table 1 Uisbs: the. known elastic: constant s 
and densities required for the- calc ul'atib n . 



TABLE I 





! Elasiic Constants 


Density 




('10^^ dynes-/cm^) 


(gm/cm-^) 




i ^11 ^iz ^ 13 - ^33: 




Zinc 


; 14*3 l-T 3:. -3: 5:: O: 4v0: 


7.1 


q 

Magnesium 


S..35 Z.53^ Z.IT 6'. 64- i;84- 


1.75 



The metals Zinc and Magnesliun arer two. different: examples 

of hexagonal material.. The dif*_fer.ence: between: and; 

in the two materials affects the response, characteristics . 

In both materials the Pure Transverse, mode. (TO has: a;- non-zero 
deviation angle for energy flow.. This: means: that: although 
the mode is a pure vibration,, associated., with. the: root' V 2 and 
lying along the v' atxis,. the energy flow is. non-zero : and not 
coincident with the x* axis (the propagation direction). 

The response for Zinc indicates that no accidental ..mode 
exists since the Quasi-Longitudinal (QL) response does not 
become zero except at the end points 0 —0° (the c direction) 
and 0 = 90° (the x-y base plane). The' absence of ' accidental 



^W. P. Mason, ed.,. Physical Acoustics v; . IIIB (New York, 

1965 ), p. 32 * 
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response can be verified by equafdon- (23). Tbee Qtaasi-- 
Transverse (QT) mode has a zero value between.. therendipoints 
but this does not imply pure raede response since^ the ; amplitude 
^ non-zero value at this point. For. Zinc: the- general 
effect of elastic anisotropy on the energy flow is.- to : channel 
the energy of most of the waves toward the base^ plane. 

An accidental pure mode is odservable in Magnesium. The 
Quasi-Longitudinal mode has a. zera value at' 6 ==51’. 6°'. 

Analysis of the amplitude shews that, both A^^- andi A^^- are 
zero along this direction, of propagation. Thus.^^ the amplitudes , 
and hence the vibrational responses, are colinearrwith the 
primed coordinate system.. The Quasi— Transverse; mode ; has 
essentially the same form of response as Zinc: but . is opposite- 
ly directed. This is due to the difference in;. the. value of 
(Cii - materials-.. The magnitude: of the 

deviation angles in Magnesium indicates that* the- material is 
essentially isotropic.. 
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Cl 

'a^ 




Oaa) ‘9 ‘'aoioaA NOiivovaoHd woaa NoiiviAaQ ao aaoNv 



24 



ANGLE OF PROPAGATION DIRECTION FROM Z-AXIS, 6, (DEG) 
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ANGLE OF PROPAGATION DIRECTION FROM Z-AXIS, 6, (DEG) 




1 

1 



TL .. cubic: symmetry; crystals 



Calaula-ttona £ar- crubxc:: ^nunetTrjr- cryBtals: follow the same 
pattern as that for heocag^onall cr.ystals'. . Ther conventional 
axes are parallel to the cr.yst-all basis.' vectors : a , ,b , and c. 
The elastic constants are =- c:|-|- .and: C^^ = 

Transformation of these const-ants: leads: to: ther equations 
shown below.. 

ll ’ 1 l h P' P' P" P' P ' P 

C^l = (K t m +- n )C^ t- m +-■ £ n +-m.- h ) (C^2 

= 2n^a~^(Jl^ +• m^ + A^)(C^^. — C'^2' 
cjg = ZlVc'^CC^- - - 2Ei|j|) +-g;,j, 

—1 4- 21 Z" (30) 

= ~na (SL + ra - n a )(C-,-,- — 0 ^^-— 2 Cj^^) 

= -Kma~^(,£.^ - m^')‘(C^ — C^- — 2C^^) 

= Jtmna“^( Jl^ - m^)'(C:^ - C'^. — 2C^|^ ) 



For propagation in the- x-^,, x-z:, . y;-s: or. (110) planes the 
constants C£g- and C£^ vanish yielding_ a, partially diagonal 
secular equation as for the hexagonal- case . The velocity 
calculations are identical to. those for. hexagonal materials. 
The relations from the secular, equations' are 



v: 



v: 



2 

2 

1,3 



= C 



66/p 



2F ^^4i ^ 



4 



(C£^ + C£^)' 



- 4 (C^C’ --c:l) ) 



'11^55 



15 



(31) 



Further response analysis' in: the, above planes yields 
amplitudes that vanish for certain modes of propagation and 
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elastic constants that vanish, hy choicBi of; pi'ane. of propaga- 
tion. Thus only C^. and: r.eraaih. tn: bef inserted : in 

the resultant energy flow equations: whihh- arer siniiiar:' to . 

Eq. (25) for hexagonal symmetry., The.- graphs: on. the; following 
pages show solutions for Aluminum with sound: velocity propa- 
gation in the (IQQ) plane (x-y plane.) and' the: (IID) plane 
(4> = ^5°)* The results agree with Muagrave:' s; data. . The-; 
elastic constants for Aluminum are 

X IQ^^ dynes'/cm^" 

C ^2 ~ ^ 10^'^ dynes/cm^- 

^44. ~ dynesycm^" 

“IT- 

The value for density is Z..6^^ gm/cm .. These: values; are : taken 
from Musgrave . 

The results for (IQQ) plane propagation, lir. Aluminum- show 
the existence of the accidental mode at 0 =- 45°" which is 
predicted by Borgnis.. The. high symmetry of: th.er resppnse 
curves is expected in the case of cubic: crystals:. Inr'addition 
the deviation angle for the Eure Transverse mode- is: zero. 

Thus the energy flow is calinear with the propagation direc- 
tion. For propagation in the (110) and (100) planes the 
angle, is zero. 

The response for (110) propagation merits, particular 
attention. The Quasi-Longitudinal mode deviation angle goes 
to zero at 6 = 54.6° which is the. [Hi]' direction. This is 
not an accidental mode but is^ indicative of:' the; fact" that 
is zero because of propagation direction. Thus the 
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secular equatlLon. rs diagonal- and. the threes modes : are pure . 

In addition, tha velacdtxesr and are: equal land: the devi- 
ation ang:les for tiieae: rorois: ar.e: equal andl nonr-zero . The 
phenomenon, of internall oonicall refraction exist s:f6r this 
propagation direction*. It. can.be shown that P^, ,and_ 

CPy.r "t ^ constant*. Since the response-' for is 
equal and op-posite to that fbr.'A2, the response for- any 
linear combination of the' two:,, (making an angle ;y> .with the 
y^-axf s) would occur at. twice: the angle y » • ih; the : -y direction. 
In other words the energy flow describes a: cone about the 
propagation, direction*. If' the: amplitude for: one-mode is 
rotated the response- rotates' in the opposite- direct ion at 
twice the angular frequency' of rotation. 



% 
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ANGLE OF PROPAGATION DIRECTION FROM Z-AXIS, 6, (DEG) 



o 
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ANGLE OF PROPAGATION DIRECTION FROM Z-AXIS, 0, (DEG) 



¥H.. TEEffAGONAL" SYMMETRY. 



Solutions for tetragannl. crystniss r.epr^sentr a- step: 
forward since there is Ix-trfc.TeL- daia-- ayalTablee to: her compared 
with. The I.R.E^ comrentxans- d;ef ihe: the: z^xls: parallellto 
the basis vector c and the x-axi's: par^lTeil to: the: basis : 
vector a, or the twofold symmetry a:xis'.. The: two: types ■ of 
tetragonal class materials are: dd'stihgui'shabler by- the^fact": 
that = 0 for Tetragonal 1 materials.:. . The: transformation, 

of elastic constants is atherwise identical.. The: transforms- - 
tions for Tetragonal 11 materials" are: 

^11 " “^^^'^11 t- 2h-^a^(C.^, +- i-nLC^^' 

+ 4ji^m^Cgg + 42,m(j,^ — ra^)C.^ 

2-2.4^ 4,-^ ^ 2T 2:- V 

= n a (£ + m + 2x. ra n ot. (C^.p +~ 2.Cgg.) 

- n^a^(2C^2 “ ^33^ i +- 4-mn'.^a"^<£^'--m^)C^g 

= 2£^m^a"^(C^^ - C^) -t n^C^ +- o:”^£.^ —nrO^Cgg, 

- 4£ma”^(£^ - in^)C-]-g (32) 

= -na“^(£^ + m^)C^^ - 2£^m^n t na(i:— 2h^MC^2 

+ n^aC^^ ~ 4£^m^na~^Cgg - 4£mna~^(£^' — ra^)C^g 
C|g = -£ma"^(£^ - m^)(C^^ - -2Cgg) 

-1. 4 ^ 4 ^„2 2.^ 

+ a (£ + m - d£ m )C-|^g 

= £mna~^(£^ - m^)(C^j^ — — ^66^ 

-2,„4 4 r2 2.^ 

-na (£ +m — 6£ra)C^g 
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Once again tlxe ctuaice oT (XQCI) ocn (HDl) pXanee propagation 



eliminates the c onatrants amd and: the: secuiarr equation, 
is partially diagonal.. Tha energy IXo'w cadhulatlDns; are: 
simplified since and ar.a identically- zero; 

The remaining constants are 
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4 s 



n » 

^35 






(l-n^)C^2 ^^^13 

n^(l-n^)CC^j^ + ( l-Sh^CX-n^-) )C-^3, 

-nd-n^)^/^ (C 3.3 - 

-n^(l-n^)^^^ - n.(l4n,^)^^^(l-ln^)C'^, 

+ n(l-n^)^'^^C33 - 2n.(X-Hi.^)^'^^(l-2h^yc^^ 
n(l-n4^''^(C„ - Cjg) 



(33) 



The resultant velocity and. energy Clov'i calculations: are: 
identical to those for hexagonal syrarae-try . . The: results: are 
shovm on the next graph for (ICfl) plane- pnopagation: in: Tin. 



The required data for calculation is given: b.elbw. . 
= 8.27 X 10^^ dynes/ cm^ 

C ^2 = 5.78 

C 13 = 3.42 ” 

C 33 '“-31 

Ciit, = 2.69 
46 = 2.82 
P =7.39 gm/cm^ 



Si 



^Ibid. 
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The response is similar fea tliat: ibr. Z3_hc: lir. that, no. 
accidental mode occurs., In addi.txorr, the. smaill deylatibn;- 
angles again indicate that the mat.erieX is: essentially 
isotropic. For the tetragonai resultsi the-: deviatlbn.. angles - 
are not reversed in sign.. 
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ANGLE OF PROPAGATION DIRECTION FROM Z-AXIS pEG) 



vnr., caNCLusiDU 



Agreement with. pubLished' data'. indidat-eB." ttiat: the method 
presented in this paper is’ a vaiid: one: and: is; particularly 
useful in the determination, ctf energy fiow ' components ; Thus 
it is possible to calculate the direction. of: energy flow for 
waves propagating in an. arbitrary direction. . The: results 
also show that, perhaps surprising^ly ,, there: may be. directions 
of propagation for which all the modes ar.e. pure but* the . trans 
verse vrave energy flow vector and: the: propagation: vector are 
not colinear. This method can.be extended: to: a: crystal ! of 
general symmetry^ The effecte ocf* crystall stability and multi 
pie mode response are areas for further. Ih'V'estigation. 
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